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further  experimental  and 

THEORETICAL  STUDIES  OF  UEDEREXPAEDED 
JETS  REAR  THE  MACH  DISC 


L.  D'Attorre  and  F.  Harshbarger 

ABSTRACT 

Results  of  experimental,  flow  visualization  studies  in  the  neighbor¬ 
hood  of  the  Mach  disc  in  underexpanded,  inviscid  jets  axe  presented.  The 
cases  include  sonic  orifices,  conical  nozzles,  Foelsch  coordinate  nozzles, 
and  some  full  scale,  in  flight  missile  data.  The  exit  to  ambient  static  * 
pressure  ratios  go  up  to  3^0  and  the  ambient  Mach  numbers  go  up  to  4.3* 

In  each  case  studied  the  method  of  Bowyer,  D'Attorre,  and  Yoshihara  appears 
to  provide  a  reasonable  technique  for  predicting  the  size  and  axial  position 
of  the  Mach  disc.  In  addition  it  is  shown  that  the  curvature  in  the  Mach 
disc  must  be  considered  if  a  consistent  triple  point  calculation  is  to  be 
made. 
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FURTHER  EXPERIMENTAL  AND 


THEORETICAL  STUDIES  OF  UNDEREXPANDED 
JETS  NEAR  THE  MACH  DISC 


L.  D'Attorre  and  F.  Harshbarger 


I.  INTRODUCTION 

Several  different  approaches  for  calculating  the  position  and  size 
of  the  Mach  disc  in  underexpanded  flows  have  been  postulated.  D'Attcire 
and  Harshbarger"'"  compared  the  results  of  a  variety  of  different  experi¬ 
mental  conditions  (including  hot  and  cold  flow  as  well  as  a  quiescent  and 

a  supersonic  ambient)  with  results  calculated  by  the  procedures  of  Bowyer, 

2  5  4  5 

D'Attorre,  and  Yoshihara,  Adamson  and  Nicholls,  and  Eastman  and  Radtke. 

1  2  R 

This  initial  comparison  indicated  that  the  Bowyer,  D'Attorre,  and  Yoshihara  ' 

approach  provided  results  in  better  agreement  with  the  varied  experimental 

k  5 

results  than  the  ether  two  approaches.  * 

Several  additional  experimental  conditions  have  been  treated  recently 
which  provide  more  justification  for  the  initial  conclusion  and  which  re¬ 
solve  several  problems  which  were  brought  to  the  attention  of  the  authors. 

It  is  the  purpose  of  the  present  report  to  document  this  additional  work. 
Another  report  provides  a  complete  summary  of  the  results  of  this  program.^ 
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II. 


SUPERSONIC  AMEIEIOT,  HOT  FLOW 


A.  Description  of  the  Experimental  Apparatus 

A  rocket  engine  with  a  Foelsch  nozzle,  an  8:1  expansion  ratio,  and 
burning  gaseous  oxygen  and  liquid  RP-1  was  studied.  The  nozzle  exit 
diameter  was  1.500  inches.  The  nominal  total  propellant  flow  rate  was 
0.695  lb/sec.  The  engine  was  operated  at  a  nominal  500  psia  chamber 
pressure  and  at  a  propellant  mixture  ratio  of  2.25  to  1  (GO^/rP-I)  by 
weight.  The  rocket  was  cooled  by  water  and  could  be  operated  for  more 
than  an  hour  without  damage.  (The  rocket  nozzle  was  designated  as  Type 
II  in  the  previous  report.^") 

The  studies  were  carried  out  in  the  supersonic  leg  of  the  AEDC  l6 '  x 
l6 '  propulsion  wind  tunnel.  The  rocket  was  fired  axially,  but  about  two 

feet  above  the  center  line.  The  products  of  combustion  were  removed  by 

a  scavanging  system.  The  rocket  jet  was  photographed  through  windows 
in  the  top  of  the  wind  tunnel  by  cameras  operated  remotely.  The  propellants 
were  not  doped  with  sodium. 

B.  Experimental  Conditions  and  Results 

The  engine  was  operated  at  a  pressure  altitude  of  50,000  ft  and  at 
an  ambient  flow  Mach  number  equal  to  1.75*  The  photograph  of  this  flow 
is  presented  in  Figure  1.  The  calculated  positions  for  the  free  stream 
shocks  the  jet  boundary,  and  the  jet  shock  structure  have  been  super¬ 
imposed  on  the  photograph  and  on  the  same  scale.  The  value  of  the  axial, 

distance  divided  by  the  nozzle  exit  diameter  (X m/dg)  to  the  Mach  disc  as 
calculated  by  the  technique  of  Bowyer,  D’Attorre,  and  Yoshihara2*  3  is  6.3. 
The  measured  distance  to  the  Mach  disc  is  6.3  ±0.15.  The  calculated 
diameter  of  the  Mach  disc  is  0.10  whereas  the  measured  diameter  of  the  Mach 
disc  is  0.1  ±  0.003. 


2 


The  flow  field  for  nearly  the  same  experimental  conditions  is  provided 

in  Figure  18.  Also  plotted  in  Figure  l8  is  the  value  of  the  pressure  behind 

the  incident  shock  as  a  function  of  the  axial  distance.  According  to  the 

5 

method  of  Eastman  and  Radtke  the  position  of  the  Mach  disc  occurs  when  this 
pressure  is  a  minimum.  The  position  so  determined  is  Xja/de  =5-9  which  is 
7$  different  from  the  experimentally  determined  value. 

The  same  rocket  engine  was  fired  also  at  a  pressure  altitude  of 
66, 500  ft  with  both  a  free  stream  Mach  number  of  2.15  and  a  quiescent  ambient. 
These  results,  along  with  theoretical  calculations  for  the  same  conditions, 
have  been  presented  previously.  However,  the  quiescent  ambient  results 
were  calculated  incorrectly  and  have  been  corrected  in  Section  III  of  this 
report. 


III. 


QUIESCENT  AMBIENT,  HOT  FLOW,  FOELSCH  (A  CORRECTION) 


In  Reference  1  consideration  is  given  to  an  experiment  in  which  the 

flow  from  an  8:1  area  ratio  nozzle  (hot  flow,  Foelsch  contour,  nozzle 

Type  II,  see  Section  II)  was  photographed  in  an  altitude  tank  at  a  pressure 

altitude  of  66,500  ft.  The  photograph  was  reproduced  as  Figure  7  is 

Reference  1  with  the  correct  dimensions  indicated.  However,  in  analyzing 

the  result  the  nozzle  exit  diameter  scale  was  misused.  In  Figure  2  of 

this  report  we  provide  ‘the  same  photograph  as  was  contained  in  Figure  7, 

Reference  1.  Figure  2  contains  lines  for  the  position  of  the  incident 

shock,  Mach  disc,  reflected  shock,  and  streamline  at  the  triple  point  as 

2  3 

calculated  by  the  method  of  Bowyer,  D'Attorre,  and  Yoshihara.  In 

Figure  lk  of  this  report  is  given  a  calculated  flow  field  along  with  a 
plot  of  the  pressure  ratio  across  the  incident  shock  as  a  function  of 
axial  position  and  the  position  of  the  Mach  disc  as  determined  by  the 

5 

experiment,  by  the  technique  of  Eastman  and  Radtke,  and  by  the  technique 

2  3 

of  Bowyer,  D'Attorre,  and  Yoshihara.  *  In  Figure  3  are  given  the  revised 
values  for  the  various  triple-point  parameters  (corresponding  to  Figure  20 
in  Reference  l). 

It  will  be  noted  in  Figure  lk  that  the  axial  position  of  the  Mach 
disc  determined  theoretically  is  X^/r^  =  21. k,  by  the  technique  of  Bowyer, 
D'Attorre,  and  Yoshihara  is  22. k,  by  the  technique  of  Eastman  and  Radtke 
is  18.2,  and  by  the  technique  of  Adamson  and  Nicholls  is  l8.k. 


QUIESCENT  AMBIENT,  HOT  FLOW,  CONICAL 
A.  Description  of  Experimental  Apparatus 

A  rocket  engine  with  a  conical  5:1  area  ratio  nozzle  burning  gaseous 

oxygen  and  sodium-doped  gasoline  was  used  to  obtain  photographs  of  an 
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underexpanded  jet  in  a  quiescent  ambient  by  Hendershot  and  Leo.  The  exit 
and  throat  diameters  were  O.366  inches  and  0.l64  inches,  respectively. 

The  nozzle  half -angle  was  15° •  The  total  propellant  flow  rate  was  0.045 
lb/sec.  The  engine  was  operated  at  a  nominal  340  psia  chamber  pressure 
and  at  a  propellant  mixture  ratio  of  1.8  (GOX  to  RP-1,  by  weight).  The 
rocket  was  uncooled  and  was  operated  for  approximately  2  seconds. 

The  studies  were  carried  out  in  a  1000  ft^  vacuum  chamber,  6-ft  in 
diameter  and  approximately  36 -ft -long.  The  rocket  was  fired  axially  and 
photographed  through  a  12- in. -diam  viewing  port.  The  altitude  tank  was 
pumped  down  to  30(41  (approximately  185,000  ft  pressure  altitude)  and  the 
engine  fired.  The  pressure  altitude  dropped  from  the  initial  value  to 
130,000  ft  during  the  first  second  of  the  rocket  run.  Stable  chamber 
pressure  was  achieved  within  0.20  to  0.25  sec  after  ignition  at  which  time 
thte  pressure  altitude  was  approximately  150,000  ft. 

A  diaphragm-type  differential  pressure  transducer  was  designed  and 
fabricated.  The  pressure  on  both  sides  of  the  diaphragm  was  equilibrated 
before  the  run  at  the  measured  pressure  altitude.  Strain  gauges  were 
bonded  to  a  cantilever  beam  which  made  contact  with  the  diaphragm.  The 
strain  gauge  output  was  calibrated  for  known,  small,  pressure  differences. 
This  pressure  transducer  was  mounted  inside  the  vacuum  chamber,  3-ft  directly 
below  the  rocket  engine.  The  strain  gauge  outputs  were  recorded  continu¬ 
ously  during  the  run. 

A  1$  solution  of  sodium  butyoxide  mixed  with  the  fuel  provided  satis¬ 
factory  self-emission  in  the  plume  for  photography  at  all  altitudes. 


A  more  complete  description  of  this  work  has  been  given  by 
7 

Hendershot  and  Leo. 

B.  Experimental  Conditions  and  Results 

Photographs  were  obtained  at  seven  different  pressure  altitudes. 
The  altitudes,  photograph  scale-factor,  and  so  forth  are  given  in 
Table  I. 


TABLE  I.  SUMMARY  OF  CONDITIONS 


Figure  No. 

Pa 

(mm  Hg) 

Pressure 

Altitude 

(feet) 

Pt/Pa 

Scale  Factor 
for  photos 
/photo  scaled 
Vmodel  scale/ 

(in. /in.) 

4 

5.50 

108, 000 

3200 

0.258 

5 

3.90 

117,000 

4510 

0.258 

6 

2-94 

125,000 

5960 

0.258 

7 

1.84 

135,000 

9560 

0.2295 

8 

1.70 

137,000 

10350 

0.254 

9 

1.60 

140,000 

11000 

0.254 

10 

1.05 

150,000 

16750 

0.278 

It  will  be 

noted  in  Figures  4  through 

10  that  the 

onset  of  the  bright 

region  indicating  the  position  of  the  normal  shock  is 

curved,  with  the 

center  of  curvature  away  from  nozzle  exit. 

In  all  schlieren  photographs 

which  have  been 

.  examined  the 

;  curvature  of  the  Mach  disc  is  such  that  the 

center  of  curvature  is  located  on  the  same 

side  of  the 

shock  as  the  nozzle 

exit.  This  difference  between  a  photograph  of  the  self-emission  and  a 

1  1 

schlieren  picture  was  noted  previously.  D'Attorre  and  Harshbarger  indi¬ 
cate  that  as  the  photographic  exposure  time  is  reduced  by  using  a  high 


6 


frame -rate  movie  camera  the  individual  exposures  show  the  shock  to  be 
flatter  than  is  observed  in  photographs  4  through  10.  Also  the  axial 
.position  of  the  shock  changes  from  frame  to  frame  in  the  experimental  condi¬ 
tions  they  studied.  It  is  believed  that  the  effect  of  a  somewhat  unsteady 
shock  position,  when  integrated  over  the  l/25  second  exposure  time  in 
Figures  4  through  10  results  in  an  apparent  inverse  curvature  or  shape  for 
the  Mach  disc.  The  selection  of  the  exact  position  for  the  stopping  shock 
in  Figures  4  through  10  was,  as  a  result,  somewhat  arbitrary. 

The  values  selected  for  the  Mach  disc  position  divided  by  the  nozzle 
exit  diameter  (X ^dg)  are  given  in  Table  II.  Also  given  in  Table  II 

SUMMARY  OF  RESULTS 

Mach  Disc  Mach  Disc 

Axial  Position  Diameter 

Vde  Vde 


4  • 

97-0 

19-0 

7-9 

5 

136.5 

24.0 

10.1 

6 

181.0 

27.0 

13.2 

7 

290.0 

31.0 

15.6 

8 

313.5 

32.0 

15.7 

9 

333.8 

37-7 

16.9 

10 

506.0 

39-5 

19.5 

TABLE  II. 


Figure  No. 


^e^e 


are  the  values  selected  for  the  Mach  disc  diameter  divided  by  the  nozzle 
exit  diameter  (d m/dg) . 

A  calculation  was  made  for  the  axial  position  of  the  Mach  disc  by 

two  methods  for  the  case  given  in  Figure  9*  In  Figure  9  is  superimposed 

the  incident  shock,  the  Mach  disc,  and  the  boundary  as  computed  using  the 

method  of  Yoshihara,  D'Attorre,  and  Bowyer  ’  and  the  pressure  behind  the 

incident  shock  as  applicable  to  the  method  of  Eastman  and  Radtke . ^  The 

I4. 

technique  of  Adamson  and  Nicholls  should  not  be  applied  to  this  case 
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as  it  lies  outside  their  bounds  of  validity.  A  summary  of  the  numerical 
results  is  given  in  Table  III.  Both  theoretical  approaches  appear  to 


TABLE  III.  COMPARISON  OF  THEORIES  WITH  EXPERIMENT 


X Jd 

nr  e 

Va, 

Experiment 

37-7 

16.9 

2 

Yoshihara,  D 'Attorre,  and  Bowyer 

40.0 

18.7 

Eastman  and  Radtke 

38.4 

19.1 

provide  reasonable  agreement;  for  this  particular  experimental  condition. 
Calculations  were  not  made  for  the  other  experimental,  conditions. 


V.  QUIESCE-NT  AMBIENT,  COLD  FLOW.  SONIC  ORIFICE 

Q 

Kawamura  conducted  a  number  of  experiments  on  underexpanded  jets 
and  attempted  to  explain  the  flow  in  the  neighborhood  of  the  triple- shock 
by  using  the  usual  triple-point  relations.  He  reported  good  agreement 
for  the  cases  in  which  the  ambient  to  total  pressure  ratio  (p  /p  )  is  high. 

He  reported  one  case  in  which  p  /p.  was  relatively  low  for  which  the  normal 

El  X 

triple-point  relations  apparently  fail.  This  case  was  brought  to  our 

2 

attention  by  Adamson.  In  order  to  determine  the  range  of  validity  of 

P  Q 

the  triple-point  approach  to  the  method  of  Bowyer,  D'Attorre,  and  Yoshihara 
this  Kawamura  case  has  been  studied  in  some  detail. 

A.  Description  of  Experiment 

Q 

Kawamura' s  experiments  were  conducted  with  a  10- centimeter-diameter 
high  pressure  chamber  which  was  terminated  with  a  flat  plate.  An  orifice 
was  placed  in  the  center  of  the  flat  plate.  The  orifice  diameter  was  22 
millimeters  and  its  circumference  was  beveled  at  45°  forming  a  sharp  edge 
facing  downstream.  The  working  fluid  was  dry  air.  It  was  exhausted  into 
the  atmosphere.  The  jet  structure  was  obtained  by  use  of  a  Hchlieren 
apparatus.  The  Kawamura  paper  studied  by  the  present  authors  did  not  con¬ 
tain  schlieren  photos  suitable  for  analysis. 

The  present  authors  attempted  to  duplicate  exactly  the  experimental 
apparatus  of  Kawamura.  Results  were  obtained.  However,  the  quality  of 
the  schlieren  photographs  obtained  were  not  as  good  as  those  given  by 
Ladenburg  and  Bershader  for  a  smaller  sonic  orifice  diameter  with  nearly 
the  same  pressure  ratio  as  the  one  which  gave  Kawamura  trouble.  The  re¬ 
sults  obtained  by  the  present  authors  and  by  Ladenburg  and  Bershader  are 

similar.  It  is  presumed  that  they  agree  also  with  the  results  obtained 

g 

by  Kawamura.  The  result  presented  by  Ladenburg  and  Bershader  was  obtained 
from  apparatus  described  by  Ladenburg,  VanVoorhis,  and  Winckler.^  The 
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working  fluid  is  dry  air.  The  sonic  orifice  diameter  is  10  mm.  Schlieren, 
shadowgraph,  and  interferometric  visualization  studies  were  made  for  pressure 

Q 

ratios  throughout  the  range  covered  by  Kawamura. 

B.  Experimental  Conditions  and  Results 

Kawamura  studied  the  triple-point  relations  for  p  /p.  equal  to  0.320, 

0.300,  0.190,  and  0.175.  It  was  at  the  lowest  ratio  of  p  /p.  (i.e.,  0.175) 

cL  Z 

that  Kawamura  found  the  triple-point  relations  invalid. 
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In  Figure  11  is  given  the  shadowgraph  result  of  Laderiburg  and  Bershader 
for  a  ratio  p  /p  =  0.175.  The  photograph  is  excellent  and  shows  clearly 
the  curvature  in  the  Mach  disc.  Also  shown  in  Figure  11  are  the  lines  re¬ 
sulting  from  a  triple-point  calculation  made  by  the  present  authors. 

The  triple-point  conditions  were  calculated  using  an  estimate  of  the 
flow  Mach  number  interpolated  from  the  flow  field  properties  given  by 
Ladenburg,  VanVoorhis,  and  Winckler.  The  Mach  number  used  is  =  3.49 
which  is  slightly  different  from  the  one  used  by  Kawamura  (M1  =  3*51).  We 
do  not  consider  this  an  important  difference  because  the  triple- shock  solution, 
in  general,  is  not  as  sensitive  to  the  value  of  the  incident  Mach  number  as 
it  is  to  the  shock  angles.  The  local  flow  angle  is  determined  from  the 
experiment  assuming  the  flow  is  source  flow.  The  flow  angle  appears  to  be 
0-^  =  8°.  We  assume  the  Mach  disc  is  normal  to  the  flow  locally  at  the 
triple  point.  The  tangent  of  the  incident  shock  with  respect  to  the  flow 
at  the  triple  shock  is  estimated  to  be  =  35°. 

The  triple  point  solution  was  obtained  using  graphs  published  in  the 
Appendix  of  Reference  1  for  y  =  1.4  and  M1  =  3*49.  The  angle  of  the  reflected 
shock  with  respect  to  the  incident  flow  angle  is  calculated  to  be  46.2°. 

In  summary,  at  the  triple  point  the  slip  line  angle  is  about  8°,  the  inclination 
of  the  incident  shock  with  respect  to  the  axis  of  the  jet  is  2rJ°,  the  angle  of 
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the  reflected  shock  to  the  axis  is  46.2°  -  (20.4-8°)  =  33.8°.  These 
angles  have  heen  traced  on  heavy  lines  in  Figure  11  and  show  excellent 
agreement  with  the  experiment,  in  contradiction  of  the  comments  of  Kawamura. 

8 

It  is  believed  that  Kawamura  failed  to  find  valid  triple-point  condi¬ 
tions  for  this  case  because  he  assumed  that  the  Mach  disc  was  everywhere 
perpendicular  to  the  axis  of  symmetry.  As  a  result  the  relative  angles 
between  the  incident, -reflected,  and  Mach  disc  shock  are  not  proper.  The 
deviation  from  perpendicularly  to  the  axis  at  the  triple  point  (about  8° 
for  this  condition)  is  sufficient  to  introduce  a  discrepancy  in  the  triple 
point  calculations.  At  higher  ratios  of  p  /p  the  Mach  disc  is  smaller, 
the  curvature  is  less,  and  the  slope  of  the  Mach  disc  at  the  triple  point 
is  much  more  nearly  normal  to  the  axis.  As  a  result,  Kawamura  was  able  to 
obtain  valid  triple-point  relations  for  these  less  underexpanded  conditions 
by  assuming  the  shock  was  everywhere  perpendicular  to  the  axis. 
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VI.  FLIGHT  DATA 


A.  Description  of  Experiment 

Excellent  tracking  photographic  coverage  is  obtained  on  missiles 
during  launch  from  the  Western  and  Eastern  Missile  Ranges  as  standard 
operating  procedure.  For  certain  camera  angles  and  missile  altitudes 
the  plume  structure  is  quite  well  defined.  The  Thor  and  Jupiter  missiles 
make  use  of  the  same  single  engine,  8:1  area  ratio,  semi -bell  nozzle. 

The  flow  angle,  static  pressure,  and  Mach  number  are  not  uniform  across 
the  nozzle  exit,  but  are  roughly  so.  Therefore,  it  makes  sense  to  study 
some  of  this  field  data.  It  provides  a  ready  source  of  jet  photographs 
with  supersonic  ambient  flow.  Unfortunately  all  the  missiles  fly  roughly 
the  same  trajectory  so  that  only  a  narrow  range  of  experimental  conditions 
(Mach  number  and  pressure  ratio)  is  obtained  by  this  technique. 

B.  Experimental  Conditions  and  Results 

A  photograph  for  one  set  of  conditions  is  presented  in  Figure  12. 

The  aspect  angle,  ambient  Mach  number,  and  pressure  ratio  used  for  the  data 
reduction  are  given  on  the  figure. 

It  will  be  noted  that  it  is  impossible  to  determine  the  exact  size 
of  the  Mach  disc  from  this  photograph.  However,  it  is  possible  to  esti¬ 
mate  an  upper  limit  for  the  diameter  of  the  Mach  disc.  The  axial  position 
of  the  Mach  disc  can,  in  most  cases,  be  located  with  fair  accuracy. 

The  values  for  the  axial  position  of  the  Mach  disc  and  for  the  upper 

limit  of  its  diameter  are  given  for  seven  cases  in  Table  TV.  In  Table 

TV  pe  is  the  nozzle  exit  static  pressure,  pa  the  ambient  static  pressure 

M„  the  ambient  Mach  number  relative  to  the  missile,  X  is  the  axial 
a  m 

distance  to  the  Mach  disc,  d  is  the  diameter  of  the  Mach  disc,  and  d 

’  m  ’  e 

is  the  diameter  of  the  nozzle  exit. 
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TABLE  IV.  SIZE  AND  POSITION  OF  MACH  DISC  TAKEN  FROM  FIELD  DATA 


Upper  limit  for 


Item 

pe/pa 

M 

a 

X  /d 
nr  e 

d  /d 
nr  e 

1 

7.80 

2.84 

7.28±0.20 

0.28 

2 

18.00 

3.55 

13.50±0.30 

0.52 

3 

19.30 

3.60 

13.65±0.30 

0.52 

4 

.  26.60 

3.79 

15.10iO.35 

— 

5 

30.90 

3.94 

15.40i0.35 

O.76 

6 

31.15 

3.98 

15.95^.35 

O.76 

7 

103.50 

4.325 

30.50i0.70 

3.1 

In  Figure  23  is  given  the  calculated  flow  field  for  nearly  the  same 
experimental  conditions  as  Item  7>  Table  IV.  Also  in  Figure  23  is  given 
the  variation  of  the  static  pressure  behind  the  incident  shock,  and  the 

position  of  the  Mach  disc  as  determined  by  the  techniques  of  Bowyer, 

2  3  5  4 

D'Attorre,  and  Yoshihara,  ’ J  Eastman  and  Radtke,  and  Adamson  and  Nicholls 

(even  though  the  latter  technique  should  not  be  applied  to  the  case  of  a 

supersonic  ambient) .  The  theoretical  values  obtained  are  3°.5>  21,0,  and 

30.0  (respectively)  for  X ^/d^.  The  experimental  value  (see  Table  IV,  Item 

7)  is  approximately  3O.5. 

k 

Surprisingly  the  technique  of  Adamson  and  Nicholls  provides  a  good 
estimate  for  the  axial  position  of  the  Mach  disc.  The  technique  of  Eastman 
and  Radtke'*  provide  a  bad  estimate.  The  method  of  Bowyer,  D'Attorre  and 
Yoshihara  provides  not  only  a  good  estimate  for  the  axial  position  but  also 
(within  the  accuracy  with  which  the  photo  can  be  interpreted)  for  the  diameter 
of  the  Mach  disc. 
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VII.  FLOW  FIELD  CALCULATIONS  AND  COMPARISONS 


The  variety  of  experimental  data  available  for  study  was  limited. 

The  results  to  date  indicate  that  the  technique  of  Bowyer,  D'Attorre, 

2  3 

and  Yoshihara  provides  a  good  estimate  for  the  size  and  position 
of  the  Mach  disc.  As  a  result,  a  number  of  flow  field  calculations  were 
undertaken  for  a  variety  of  conditions  for  the  purpose  of  comparing  with 
the  other  two  theoretical  techniques^ 

Theoretical  calculations  were  made  for  the  inviscid  flow  field 
structure  using  the  Bowyer  program. ^  The  three  techniques  for  esti¬ 
mating  the  axial  position  (and  size,  where  applicable)  of  the  Mach  disc 
were  applied  to  these  representative  cases.  The  calculations  were  made 
assuming  the  flow  is  everywhere  parallel  to  the  axis  at  the  nozzle  exit, 
that  the  nozzle  exit  Mach  number  is  equal  to  3.185,  and  that  the  jet 
V  =  1.225  and  the  ambient  y  =  1.4. 

In  Figures  13  through  23  are  provided  the  calculated  lines  of  constant 

Mach  number  and  the  position  of  the  various  shocks  and  boundaries  for 

several  different  arbitrary  operating  conditions.  Also  given  is  the 

variation  of  the  pressure  ratio  across  the  incident  shock,  the  position 

of  the  Mach  disc  and  its  size  as  calculated  by  the  method  of  Bowyer, 

2  3 

D'Attorre,  and  Yoshihara,  ’  the  axial  position  of  the  Mach  disc  as  com¬ 
puted  by  the  method  of  Eastman  and  Radtke, ^  and,  where  calculated,  the 
axial  position  of  the  Mach  disc  as  calculated  by  the  method  of  Adamson  and 
Nicholls.^ 

An  analysis  of  Figures  13  through  23  indicates  that  the  technique 

4 

of  Adamson  and  Nicholls  provides  values  for  the  axial  position  of  the 
Mach  disc  wh  oh  are  in  agreement  with  those  predicted  by  the  technique  of 
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Bowyer,  D'Attorre,  and  Yo shihara  up  to  total  to  ambient  static  pressures 

(P4./p„)  equal  to  650  for  a  quiescent  ambient,  but  at  greater  pressure 

ratios  the  two  techniques  do  not  necessarily  provide  the  same  answer. 

The  techniques  of  Eastman  and  Radtke^  and  Bowyer,  D'Attorre,  and 
2  5 

Yoshihara  ’  appear  to  give  similar  results  for  the  axial  position  of 
the  Mach  disc  for  a  static  ambient  but  quite  different  results,  in  some 
cases,  when  there  is  a  supersonic  ambient. 
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VIII.  CONCLUSIONS 


The  importance  of  considering  the  curvature  of  the  Mach  disc  in 
triple-point  calculations  has  been  indicated  in  analyzing  the  experi¬ 
mental  results  of  Kawamura  for  fairly  large  exit  to  ambient  pressure 
ratios  (pe/pa) • 

Additional  experimental  evidence  has  been  obtained  which  indicates 
that  the  method  of  Eastman  and  Radtke  cannot  be  relied  upon  to  provide 
a  correct  answer  for  the  size  and  axial,  position  of  Mach  disc  for  all 
operating  conditions. 

No  experimental  conditions  have  been  analyzed  in  which  the  technique 
of  Bowyer,  D'Attorre,  and  Yoshihara  does  not  provide  a  reasonable  estimate 
for  the  size  and  axial  position  of  the  Mach  disc. 
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FIGURE  1.  Photograph  of  a  Rocket  Operated  with  Pe/pa  =  4.68  and  a  Freestream  Mach 
Number  =  1.75*  Superimposed  are  lines  corresponding  to  the  calculated 
values  for  the  ambient  shock  (white  lines),  jet  boundary,  jet  shock, 
and  Manh  disc  fblack  lines) . 
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FIGURE  3.  Triple-point  Results  (Revised)  for  a  Rocket  Operated 
with  pe/pa  =  IO.3  and  a  Static  Ambient  (see  the  photo¬ 
graph  in  Figure  2).  This  case  was  presented  incorrectly 
in  Reference  1. 
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FIGURE  4.  Photograph  of  Flow  from  a  15°  Half -Angle  Conical  Nozzle 
in  a  Static  Ambient  with  pe/pa  =  57.0. 


FIGURE  5*  Photograph  of  Flow  from  a  15  Half -Angle  Conical  Nozzle 
in  a  Static  Ambient  with  p  /p  =  136.5. 
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FIGURE  6.  Photograph  of  Flow  from  a  15°  Half -Angle  Conical  Nozzle 
in  a  Static  Ambient  with  p  ?p  =  l8l.O 


FIGURE  8.  Photograph  of  Flow  from  a  15  Half -Angle  Conical  Nozzle 
in  a  Static  Ambient  with  p  /p  =  313.5- 
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figure  11. 


Shadowgraph  of  an  Axially  Symmetric  Air  Jet  Flowing  from  a  Round  Orifice  of  10  mm 
Diameter  at  a  Reservoir  Gauge  Pressure  of  TO  lb/ in.  (reproduced  from  Reference  9)* 
Superimposed  are  calculated  slopes  for  the  incident  shock,  reflected  shock,  Mach 
disc,  and  slip  line  locally  at  the  triple-point  (by  the  technique  of  Bowyer, 
D'Attorre,  and  Yoshihara) . 
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